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Majorana Fermions, strange particles that are their own antiparticles, were predicted in 1937 and 
have been sought after ever since. In condensed matter they are predicted to exist as vortex core 
or edge excitations in certain exotic superconductors. These are topological superconductors whose 
order parameter phase winds non-trivially in momentum space. In recent years, a new and promising 
route for realizing topological superconductors has opened due to advances in the field of topological 
insulators. Current proposals are based on semiconductor heterostructures, where spin-orbit coupled 
bands are split by a band gap or Zeeman field and superconductivity is induced by proximity to 
a conventional superconductor. Topological superconductivity is obtained in the interface layer. 
The proposed heterostructures typically include two or three layers of different materials. In the 
current work we propose a device based on materials with inherent spin-orbit coupling and an 
intrinsic tendency for superconductivity, eliminating the need for a separate superconducting layer. 
We study a lattice model that includes spin-orbit coupling as well as on-site and nearest neighbor 
interaction. Within this model we show that topological superconductivity is possible in certain 
regions of parameter space. These regions of non-trivial topology can be understood as a nodeless 
superconductor with d-wave symmetry which, due to the spin-orbit coupling, acquires an extra 
phase twist of 2n. 



Back in 1937, Ettore Majorana found particles that 
arise as real solutions to the Dirac equation. These solu- 
tions, called Majorana fermions, partially obey fermionic 
statistics. While different Majorana fermions anti- 
commute, each Majorana fermion is its own anti particle. 
The creation operator of a Majorana fermion is also its 
annihilation operator or, in other words it is an equal 
superposition of regular fermionic creation and annihila- 
tion operators. This exceptional property captured the 
imagination of many and the quest to find the Majorana 
fermion began. 

In the context of high energy physics it is speculated 
that the neutrino might in fact be a Majorana fermion. 
The testing of this claim, which originated from Majo- 
rana himself, has not been possible in the past, as it 
requires a large collider like the LHC and is still an open 
question at the time of writing. 

Regardless of the nature of the neutrino and possible 
other elementary particles of the Majorana type, Majo- 
rana fermions may be realized in condensed matter sys- 
tems. In condensed matter, excitations are not limited 
to elementary particles since they may be emergent par- 
ticles that are dressed by the medium and interactions in 
the many body ground state. As such, it is conceivable 
that emergent excitations be their own anti-particles. 
Furthermore, in condensed matter anti-particles are pro- 
vided by holes in energy bands and the superposition of 
particles and holes is possible. Such a superposition of 
particles and holes occurs as an excitation in any super- 
conductor, since the number of particles is not conserved 
due to the presence of a pair condensate. 

In order to realize Majorana fermions, a system should 
exhibit pairing between two particles of the same spin. 
This requires triplet pairing and particularly a complex p- 
wave order parameter is desirable. It has been shown [T 



that topological, spin-triplet, Px + Wy superconductors 
will support Majorana fermions in their vortex cores [2]. 
Some materials have been found to have triplet p-wave 
pairing; however their topology has yet to be proven to be 
non-trivial^. Therefore, current efforts to realize Majo- 
rana fermions have had to focus on devices which lead to 
quantum states that are either topological superconduc- 
tors or analogous to them.[Tl|4H3 An interesting analogue 
of a topological superconductor was proposed to describe 
the fractional quantum Hall state in some fractions P and 
some progress in this direction has been made[8'-n!0l. In 
that state, Chern-Simons dressed particles minimize their 
interaction energy by creating a condensate whose sym- 
metry is Px^Wy Proving beyond doubt the existence of 
this state as well as the detection of Majorana fermions 
therein still remains a challenge. 

Recently, inspired by advances in topological insula- 
tors, another route to topological superconductivity has 
opened. Fu and Kane [4] have shown that a three di- 
mensional topological insulator layer placed in proximity 
to a conventional 5-wave superconductor develops topo- 
logical superconductivity. The pairing in the system is 
induced by proximity effect while the topology is inher- 
ited from the topological insulator. This occurs since the 
pairing function is projected to one of the spin-orbit cou- 
pled bands. In order to accommodate the unique spin 
structure (and Chern number) of the topological insula- 
tor, the induced order parameter must wind its phase by 
2tt in momentum space. 

The idea of Fu and Kane was further developed by 
Sau et a/. [5] in the context of semiconductor heterostruc- 
tures. Sau and coworkers envisioned a semiconductor 
quantum well with intrinsic Rashba spin-orbit coupling 
(SOC) where the Fermi surface lies in the band. The 
required gap between the two spin-orbit coupled bands 
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is provided by an out of plane Zeeman field of an at- 
tached ferromagnetic (FM) insulator layer. Meanwhile, 
superconductivity is induced by proximity to a supercon- 
ducting layer attached to the other side of the quantum 
well. 

Recently, Alicea[6 explored the possibility of elimi- 
nating the ferromagnetic insulator layer of the Sau et 
al model. Instead of the ferromagnetic insulating layer, 
Alicea suggested using a quantum well with both Dressel- 
haus and Rashba SOC while applying an in plane mag- 
netic field. The Dresselhaus SOC tilts the plane in which 
the electron spins tend to align so that the applied mag- 
netic field can open a gap, eliminating the need for the 
FM insulator and thereby reducing the complexity of the 
device. 

Another suggestion for the realization of Major ana 
fermions has been made by Cook and Franz [7 . Their 
work involves inducing superconductivity in nanotubes 
and nano-ribbons; the relation between the type of re- 
sulting order parameter and symmetries of this model 
was explored by Fu and Berg [11]. 

The key aspects of realizing an effective p-\-ip state in 
these previous devices has been the proper combination 
of SOC, band gap or Zeeman splitting and proximity in- 
duced pairing. The spin orbit coupling is responsible for 
the non-trivial spin texture, whereas the Zeeman field (or 
an intrinsic mass term) splits the bands such that only 
one of them is relevant at low energy. Meanwhile, the 
superconductor responsible for inducing pairing through 
proximity is of a simple singlet type. 

In this Letter we use these key ingredients to address 
the question of whether a topological superconductor can 
be achieved without proximity effect. In place of proxim- 
ity induced pairing we consider pairing driven by inter- 
actions. Using a variational mean field approach, both 
phases of trivial and topological superconductivity are 
found in the model studied. The topological state we find 
can be described by a superconductor with a Gtt phase 
winding which is a result of the / = 2 d-wave phase wind- 
ing and a p + projection function. 

In order to test whether interactions may lead to super- 
conductivity in spin-orbit coupled materials we consider 
a two dimensional square lattice model. The Hamilto- 
nian of the system reads 

H = Hke + Hso + V, (1) 

where the kinetic energy term Hxe is given by hopping 
on nearest neighbors. 

Hke = E (4cjv + c]^Ci^). (2) 

Here t is the hopping amplitude, are nearest- 

neighbor lattice sites and a is a spin index. The spin- 
orbit coupling part of the Hamiltonian is given by 

i^so = ^V^lHkV^k, Hk = a-dk (3) 

k 



where ^/^k = (ckt? <^ki)^, cr is a 3-vector of Pauli matrices 
and dk = {As'mkx, Asmky,2B{coskx + cos ky — 2) -h 
M) with A, B and M material parameters. This term 
can be viewed as the lattice version of the continuum 
model introduced by previous authors [4]-[6] while its form 
resembles one of the sectors of the model introduced by 
Bernevig et al. to describe HgTe quantum wells [12 . The 
finite Zeeman splitting M may be the result of a band 
gap, an external magnetic field or a magnetic field of a 
near by ferromagnetic layer. 

We choose to model the interactions with effective on- 
site repulsion and nearest neighbor attraction, such as in 
the extended Hubbard model given by 

^ = ^0 ^ rii^riii + ^ ni^rija' , (4) 

with Uq > (repulsion) and Vq < (attraction). The 
motivation behind introducing an attractive Vb stems 
from studies of a similar model without spin-orbit cou- 
pling in the context of the cuprates.[13l |14j In the case 
of the cuprates, it has been shown that on-site repul- 
sion treated in the Eliashberg formalism leads to d-wave 
pairing on bonds. This same result can be obtained in 
mean field if one introduces an effective nearest-neighbor 
attraction as in our model above. 

In order to map the phase diagram of the model in 
Eq. [l] we adopt a variational mean-field theory. Our 
method involves obtaining a variational wave function 
that is a solution to an auxiliary quadratic Hamiltonian. 
This auxiliary Hamiltonian contains the kinetic and spin 
orbit coupling parts of the Hamiltonian in Eq. [l] as well 
as a variety of order parameters. These order parame- 
ters represent all possible mean field states such as den- 
sity waves, superconductivity etc. The mean field ground 
state is found by minimizing the expectation value of the 
interacting Hamiltonian {i.e. Eq. [T]) with respect to the 
parameters of the variational wave function. These pa- 
rameters are essentially the magnitudes of the various 
order parameters of the model. The advantage of this 
method over the usual self-consistent mean-field theory 
is that it does not assume a priori the dominance of any 
order parameter. More on its application can be found 
in Reference 

To this end, we use the following auxiliary Hamiltonian 
^Aux = ^^^3,Ak^k (5) 

k 

where we have defined the 8-spinor ^k = 

^ t t t t \T 

(here Q = (7r,7r)) and the matrix 

/ h{k) A(k) \ 
Ak = . , (6) 

VA(kr -H-ky 
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FIG. 1: Sample plot of order parameters. In this figure the 
magnitude of the order parameters is in units of t and we have 
fixed A = 0.25t, B = .5t, M = At and Uo = 2t. Circles (blue 
online) are 

Ad), 

squares (black online) are A^^\ diamonds 
(orange online) are A*^^-*, x's (red online), A*^^-* and triangles 
(brown online) represent S. This simulation was done on a 
100 X 100 square lattice. The graph shows the development 
of d-\- id order since both A^^-* and A^^-* become non-zero at 
the critical coupling. 



represents the Nambu space of particles and holes. Its 
entries are 4x4 matrices: 



hik) 

and 



-Sffz 



-Sa, ?^(k + Q) 



, ?i(k)=ek + ?^(k), (7) 



A(k) = iA]f^ay i 



' Tz- 



(8) 



where = — 2t(cos + cos ky) is the tight binding spec- 
trum and Hk is as defined in Equation ([3|. In our no- 
tation ai are Pauli matrices that act on the spin and 
are Pauli matrices which denote the sublattice degree of 
freedom. 

In the above auxiliary Hamiltonian we have allowed 
for the possibility of antiferromagnetism (AF) through 
the Neel order parameter 5, as well as several channels 
of superconductivity through the order parameter Ak = 
A'^-^^ (cos kx — cos ky) + zA*^^^ sin k^ sin ky + A*^^^ (cos k^ + 
cos ky) + A*^^^. We can now find the variational en- 
ergy numerically for a given set of order parameters and 
then minimize with respect to these parameters. This 
amounts to finding the mean field ground state energy 
and wave function of the system. A representative plot 
of these order parameters appears in Fig. [l] 

As A^^) = A^"^) = throughout the plot in Fig. [ij we 
conclude that this region of parameter space does not 
support 5-wave nor extended 5-wave superconductivity. 
We also see that 5 = in Fig. 1 and so AF is also 
not a dominant order in this region of parameter space. 



The lack of s and extended 5- wave superconductivity is a 
general characteristic of the phase diagram of this model, 
however depending on how we tune the SOC parameters 
it is possible to find a state where AF (and not super- 
conductivity) is the dominant ground state. 

To demonstrate the differing ground states of our 
model we present two separate slices of the phase dia- 
gram. Figure 2 gives two plots in a space of the interac- 
tion parameters Vq and Uq; one slice is for values of A, B 
and M where the ground state is found to be supercon- 
ductivity, while the other has the SOC parameters tuned 
so that we see a phase with an AF ground state. The 
d-\-id phase in Fig. [2] is the most interesting for our pur- 
poses as it is fully gapped and therefore its topological 
invariant is well defined. It is obtained when both A^^^ 
and A^^^ are non-zero. We view this state as having a 
projected superconducting order parameter whose phase 
winds by Gtt in momentum space. An of the winding is 
due to its (i-wave nature and the remaining 27r are the 
result of the projection on one of the spin-orbit coupled 
band. We focus to this d-\-id region of the phase diagram 
and investigate the topology of this phase. 

Before ending this section, we make a few remarks in 
passing. First, Fig. [2] showcases a phase that is purely 
d-wave in nature. One might ask why we are not inter- 
ested in the topology of this phase; the reason is that 
this d-wave order parameter has nodes and as a result 
it is difficult to properly define a topological invariant 
in this case. Nevertheless, the topology of this supercon- 
ductor may be an interesting topic for further studies and 
could be found to be non-trivial. On the other hand, due 
to nodal excitations any Majorana fermions that may be 
produced will not be protected against scattering. Sec- 
ond, in the interest of preforming an exhaustive search 
for competing ground state order in this model, we have 
checked that neither spin nor charge density waves give 
a dominant ground state contribution to the model pre- 
sented here in the studied parameter regime. 

To study the topology of the d -\- id region we calcu- 
late the TKNN number [16 (equivalent to the first Chern 
number) using our optimized mean- field wave function. 
This involves selecting a region in the d ^ id phase in 
Fig. [2] and then calculating [17] 



1 

2^ 



d^kT{\i) 



(9) 



where the Berry curvature, J^, is defined using the eigen- 
states Ak|0n(k)) = £;n(k)|0n(k)) ^418 



n m^n 



\^\4>m){4>m\'^\4>n) 
{En — Em)'^ 



(10) 

where, for the sake of brevity, we have dropped the func- 
tional dependence on k, the primed sum is a sum over 
filled bands, the e*-^ tensor has the values e^'^ = — e^'^ = 1 
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FIG. 2: Slices of the phase diagram. The left plot (plot (a)) showcases a region of parameter space where the ground state of 
the model is superconductivity. We have set A = 0.25t, B = — 0.45t and M = — O.lt in this plot. The regions of the phase 
diagram are labeled A (teal online) for normal (or non-superconducting state), B (grey online) for a d-wave ground state and 
C (beige online) for d -\- id superconductivity. The right plot (plot (b)) shows a region of the phase diagram where we can see 
either antiferromagnetism or superconductivity. Here we have sets A = 0.25t, B = — 0.9t and M = — 0.05t and our labeling is 
D (teal online) for the normal state, E (grey online) for a d-wave ground state, F (green online) for d -\- id superconductivity 
and G (blue online) for AF ground state. 



and e*'* = and summation over the repeated indices i 
and j is implied. 

By calculating this invariant we can classify the topol- 
ogy of the d-\-id region as either trivial (regions for which 
we find / = 0) or non-trivial (regions for which / = 1). 
Our results are summarized in Figure 3. As the topology 
of the system is intimately related to the number of Fermi 
surfaces before interactions are turned on[TT], Fig. 3 also 
shows a sample of the Fermi surface in each topological 
region. 

Figure [3] displays one of our main results; our model 
has regions of d -\- id topological superconductivity. As 
an example, Fig. |3] shows that the d -\- id state in Fig. 
2a is topologically trivial while the like in Fig. 2b is 
topological. Further, we see from the figure some de- 
vices might exist in the topologically non-trivial region 
at some set value of M. It may then be possible to move 
the system into a topological phase by changing M via 
either an applied field or proximity to a magnetic layer. 
In this way, our results suggest that properly applying a 
Zeeman field {i.e. tuning M) to spin-orbit coupled su- 
perconductors may result in the transition of an ordinary 
superconductor to a topological one. 

In a real material, it is not possible to tune the interac- 
tion. However, this is routinely done in cold atoms. It has 
been demonstrated recently that spin-orbit coupling may 
be simulated in cold atoms [19 . This may lead the way 
for simulating topological insulators [2Ql-[22] . Our work 
suggests that if the spin-orbit coupling and the interac- 
tions are tuned correctly, a topological superconductor 
may be simulated as well. 

In summary, we have proposed a model of interacting, 
spin-orbit coupled, Zeeman split electrons on a square 
lattice. We have shown that in regions of parameter 
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FIG. 3: Topology of the d-\-id ground state phase. In this fig- 
ure we have set A = 0.25t. The topologically trivial phase is 
labeled T (beige online) while the non-trivial phase is labeled 
N (teal online). The insets show an example of the Fermi 
surface (in the first Brillouin zone) for each phase before in- 
teractions are turned on. The lower inset, corresponding to 
the d -\- id state of Fig. 2b, shows a single Fermi surface for 
the topologically non-trivial phase while the upper inset (the 
d + id state in Fig. 2a) shows two Fermi surfaces in the topo- 
logically trivial region. 



space the ground state of the proposed model is either 
di d -\- id superconductor or an antiferromagnet. Nar- 
rowing our focus to the d + id, region we have shown 
that our system supports phases of non-trivial topology. 
The topological regions in our phase diagrams exhibit 
superconductivity that exhibits Gtt winding of its order 
parameter phase in the Brillouin zone. In the same way 
that di p -\- ip superconductor may support the existence 
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of Majorana fermions in vortex cores or on edges, this su- 
perconductor will support their existence. As our model 
considers superconductivity driven by interactions rather 
than the proximity effect, it may serve as a possible sim- 
plification for device design. Finally, our work supplies 
strong evidence that Majorana fermions might be real- 
ized in certain spin-orbit coupled superconductors under 
the proper application of a Zeeman field. 

The authors would like to thank J. Alicea, E. Berg, B. 
A. Bernevig, A. A. Clerk, M. Franz, G. Gervais and S. 
Sachdev for useful conversations. This work was sup- 
ported by Natural Sciences and Engineering Research 
Council of Canada (AF, TPB) as weh as the McGih Tom- 
linson Fellowship program (AF). 



[1] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
[2] V. Gurarie and L. Radzihovsky, Phys. Rev. B 75, 212509 
(2007). 

[3] C. Kallin and A. J. Berlinsky, J. Phys.: Condens. Matter 

21, 164210 (2009). 
[4] L. Fu and C. L. Kane, Phys. Rev. Lett. 096407 (2008). 
[5] J. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 

Phys. Rev. Lett. 104, 040502 (2010). 
[6] J. Alicea, Phys. Rev. B 81 (2010). 

[7] A. Cook and M. Franz, Phys. Rev. B 84, 201105 (2011). 



[8] P. Bonderson, A. Kitaev, and K. Shtengel, Phys. Rev. 

Lett. 96, 016803 (2006). 
[9] A. Stern and B. L Halperin, Phys. Rev. Lett. 96, 016802 

(2006). 

[10] M. Dolev, M. Heiblum, V. Umansky, A. Stern, and 

D. Mahalu, Nature 452, 829 (2008). 
[11] L. Fu and E. Berg, Phys. Rev. Lett. 105, 097001 (2010). 
[12] B. A. Bernevig, T. L. Hughes, and S. C. Zhang, Science 

314, 1757 (2006). 
[13] D. J. Scalapino, Physics Reports 250, 329 (1995). 
[14] S. Onari, R. Arita, K. Kuroki, and H. Aoki, Phys. Rev. 

B 73, 014526 (2006). 
[15] T. Pereg-Barnea and G. Refael, Phys. Rev. B 85, 075127 

(2012). 

[16] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 
M. den Nijs, Phys. Rev. Lett. 49, 405 (1982). 

[17] M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. B 
82, 134521 (2010). 

[18] M. Berry Proc. R. Soc. Lond. A 392, 45 (1984). 

[19] Y.-J. Lin, K. Jimnez-Garca, and L B. Spielman, Nature 
471, 83 (2011). 

[20] F. Mei, S.-L. Zhu, Z.-M. Zhang, C. H. Oh, and N. Gold- 
man, Phys. Rev. A 85, 013638 (2012). 

[21] N. Goldman, L Satija, P. Nikolic, A. Bermudez, M. A. 
Martin-Delgado, M. Lewenstein, and L B. Spielman, 
Phys. Rev. Lett. 105, 255302 (2010). 

[22] A. M. Essin and V. Gurarie, Phys. Rev. B 85, 195116 
(2012). 



